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Abstract 

Scalar kinks propagating along the bulk in warped spacetimes provide a thick brane realisation of 
the braneworld. We consider here, a class of such exact solutions of the full Einstein-scalar system 
with a sine-Gordon potential and a negative cosmological constant. In the background of the kink 
and the corresponding warped geometry, we discuss the issue of localisation of spin half fermions 
(with emphasis on massive ones) on the brane in the presence of different types of kink-fermion 
Yukawa couplings. We analyse the possibility of quasi-bound states for large values of the Yukawa 
coupling parameter (with v, the warp factor parameter kept fixed) using appropriate, recently 
developed, approximation methods. In particular, the spectrum of the low-lying states and their 
lifetimes are obtained, with the latter being exponentially enhanced for large u^f- Our results 
indicate quantitatively, within this model, that it is possible to tune the nature of warping and 
the strength and form of the Yukawa interaction to obtain trapped massive fermion states on the 
brane, which, however, do have a finite (but very small) probability of escaping into the bulk. 
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I. INTRODUCTION 

Extra dimensions a la Kaluza-Klein 111 M have been around as an elesrant theoretical 

n 

construct for decades. Their inevitable presence in superstring theories [3||, is by now, well- 
established. Of late, warped (nonfactorisable) spacetimes have been in vogue ever since 
■landall-Sundrum introduced them in the context of the now-popular braneworld models 



a 



In such models, extra dimensions exist, can be compact or non-compact, but the usual 
schemes of 'compactification' are abandoned. The four dimensional world (the so-called 
3-brane) is visualised as an embedded hypersurface in the five dimensional bulk. Effective 
field theories on the 3-brane are obtained without invoking 'compactification' in the usual 
sense of Kaluza-Klein. Further details of how this is achieved can be found in ^ or in 
subsequent review articles jj]. Much work has been carried out on diverse aspects of such 
models in the last four years. These include the attempted resolution of the heirarchy 
problem Q a, questions about the loeahsafo,. of v^ous types of fields on the brane L 
particle phenomeonology in the braneworld context [J] and cosmological consequences 
The actual existence of warped extra dimensions as well as a firm foundational basis for these 
models however, still remain open issues. Phenomenological work supported by experimental 
verification and simultaneous progress in formal theory (via superstrings or any other yet-to- 
be-discovered fundamental unified theory) will perhaps be able to provide conclusive answers 
about the future of the braneworld scenario. 

In the early braneworld models, the bulk spacetime was assumed to have only a negative 
cosmological constant (anti-de Sitter space in five dimensions). Subsequently, we have had 

models with fields in the bulk-the simplest of them being those with a bulk scalar field 

□ 

dependent only on the fifth (extra) coordinate j9[. In this paper, we first consider exact 
solutions of the Einstein-scalar equations with a sine-Gordon potential. There exists a full 
solution of this system with the scalar field configuration being a kink. Such a configuration, 
the so-called thick brane, provides a realisation of the braneworld as a domain wall in the 
bulk, a fact which has been illustrated in several examples in the literature Q]. The warped 
background spacetime, in which the bulk sine-Gordon soliton is a self-consistent solution, 
turns out to have a non-constant but asymptotically negative Ricci curvature. 

In the scenario mentioned above, where we have an exact solution of the full Einstein- 
scalar system, we address the question of localisation of fermions with emphasis on the 
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massive ones. Though largely a toy model, we are able to find some interesting analytical 
results, which add to the only work on massive fermions in the literature To address 
the issue of localisation, we couple the fermion field with the scalar field through a Yukawa 
coupling, as has been done in many articles on fermion localisation. Subsequently, we 
investigate and compare the localisation effects for two different types of Yukawa couplings 
in the background of this scalar sine-Gordon kink in the warped line element. Known 
results for the massless fermions are verified. Using approximate methods, we are able to 
delineate which Yukawa coupling can yield better results (i.e. have more trapped states). In 
addition, we make use of a new, recently developed formalism for estimating the lifetime of 
the states and show how the lifetime can indeed be exponentiallylarge for large values of the 
Yukawa coupling parameter and moderate values of the warp factor parameter. In summary, 
the available exact solution for the scalar field and the warp factor helps us to understand 
analytically, the nature of the effective potential, the trapped states, their lifetime and the 
comparative role of different Yukawa couplings. 



II. THE EXACT SOLUTION FOR THE WARP FACTOR AND THE SG KINK 

We begin by writing down the Lagrangian and the equations of motion for Einstein gravity 
with a cosmological constant, minimally coupled to a real scalar field with a potential V{(j)). 
The action for such a system is given as : 



S 



{R-2K)--g^'da(t>d,<j)-V{<P) 



-gd^x (2.1) 



where gab is the five dimensional metric tensor with signature (- + + + +) and R is the 
Ricci scalar. A is the five dimensional cosmological constant, = SttG^ where G5 is the 
five dimensional Newton constant. 

We assume a line element, in the five dimensional bulk, of the form : 

ds^ = da^ + e-^f^^^ijdx'dx^ (2.2) 

where a is the fifth coordinate and /(cr) is the warp factor. The scalar field is considered 
to be a function of a only. With these ansatze the Einstein scalar system reduces to the 
following system of coupled, nonlinear ordinary differential equations : 
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f" = a<p" (2.3) 



r = ^ (r - - 6 (2-4) 

(2.5) 

where a = k^/S. The first two of the above set are the Einstein equations and the third 
is the scalar field (Klein-Gordon) equation, which follows from the Einstein equations and 
is not independent. 

The obvious question is : is there an exact analytical solution to the system of Eqns (2.3)- 
(2.5)? For the Higgs potential, Ringeval et al have numerically investigated the solutions 
and found kinks and the corresponding warp factors. We shall now show that the Eqns 
(2.3)-(2.5) can indeed be solved exactly for a sine-Gordon potential {V{(j)) = B (^1 + cos ^ . 
Variants of this solution with different choices of V{(f)) have appeared in the literature (l2|. 
We comment on this at the end of this section. We now write down the solution explicitly 
below. 

m = ^Mlncosh(^a) (2.6) 

0(a) = 2Atan-^ fexp —a) - — (2.7) 
V a J 2 

The constant B in the sine-Gordon potential is given as : 
B = ^{a + ^).Also., = ^^. 

It is important to note the presence of a negative cosmological constant in both the warp 
factor and the soliton. In addition, we emphasize that the kink/soliton is constructed in 
such a way that 0(cr) is an odd function about the origin. The second term in the expression 
for 0(cr) is required to achieve this. 

This model with a bulk SG potential provides a 'thick brane' realisation of the Randall- 
Sundrum scenario where the SG field and its soliton configuration dynamically generate this 
domain wall configuration in the background warped geometry. In addition, as is obvious 
from the functional from of /, there is no discontinuity in the derivative of / at the location 
of the brane. The warp factor is smooth everywhere and has all the necessary features. 

The line element in the warped bulk space time is given by 



ds^ = da^ + cosh"^'' {—a)r]ijdx'dx^ (2.8) 



where z/ = —\r-^ 

Kl V D 



Notice that the metric is completely non-singular for the full domain of the fifth coordi- 
nate. It describes a space of negative Ricci curvature given by 



20 ) tanh^ —a 



(2.9) 



_av42 VaA2 J a 
It is straightforward to check that the above function for R is not always negative with 
a asymptotic (a ±oo) value of — ^|A|. One cannot get constant negative Ricci curvature 
for any real value of the parameter A. In the above solution, the parameter A is free whereas 
B is fixed in terms of A, |A| and a. 

Furthermore, since the above f{a) and 0(cr) consistently solves the Einstein-scalar system 
we may want to know precisely the energy density and pressures which generates such a line 
element. These are given below. 

'^./^sech^^a-IAItanh^^. 
a V 6 a a 



'Px,y,z = 3— — sech —a — |A| tanh —a (2-10) 



p,^ = lAltanh^— cr (2.11) 
a 

Note that, asymptotically, the stress energy is that of anti-deSitter space in five dimen- 
sions. At cr = (i.e. the location of the brane) we have Px,y,z = —p and = - an effective 
cosmological constant on the 3-brane. The energy density and pressures are plotted in Fig- 
ure 2. It is worth noting that the matter stress-energy which acts as a source for the warped 
geometry satisfies the Null Energy Condition (NEC ^ p + Pi > 0) though the Weak and 



Strong Energy Conditions (WEC and SEC) are violated (WEC ^ p > 0, p + Pi > 0) 

It is worth noting that there exists related work on similar solutions for (p and /. In a 
most general case one might consider the following : 



f{a) = Ailn(cosh6a) ; 0((t) = 2A2tan-^(exp(6(T)) - — ^ ; 1/(0) = C + Dcos(2^) 

2 A2 

(2.12) 

Using these in the equations of motion one obtains the following constraints on Ai, A2, 
C, D and |A|. 
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FIG. 1: Variation of the energy density p (gray line) and pressure ( pi (black line) and Pa (dark 
gray line)) with a corresponding to the warped geometry given by the metric (4.3) 



Ai = aAi ; Aib' [1 + 4A(j = 8B (2.13) 
C = -^ + ^ + U (2.14) 

+ (2.15) 

It is easy to note from the above that if A = it is not possible to have C = D. Various 



expressions for V{(f)) with the corresponding solutions have been noted in the literature 
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though these are all for A = 0. The solutions for |A| = and Ai ^ A2 can be matched with 
the solutions for |A| 7^ with Ai = A2 = 2D, if, in the former, we assume A2 — Ai = 2^-^. 
This means that an 'effective' A will appear within the scalar potential. However, we prefer 
to keep A separate. As we shall see later, the mass spectrum of the localised fermions will be 
dependent onrip and | A| . It is also worth mentioning that in some of the literature on similar 
domain wall solutions the SG soliton is taken to be a neither-odd nor-even function. We 
have used a shifted version (shift along the axis, in order to maintain the odd character of 
the soliton solution and also to have = at the brane location). The use of the purely odd 
soliton is important because, as we note in the next section, the symmetry of the effective 
potential for massive fermions w.r.t. the brane location is respected only if (p is purely 
odd. This is desirable because the physics on either side of the brane should have identical 
features {Z2 symmetry). 
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In this context, one must also mention the extensive work done in the recent past on 
non-supersymmetric as well as supersymmetric , domain walls in diverse models and 
theories. 

Let us now focus on the question of localisation of massive and massless fermion fields 
on the 3-brane. 

III. FERMION LOCALISATION 

As with other fields such as scalars, spin one fields, gravitons, gravitinos, it is a pertinent 
question to ask whether spin-half fermions are/can be localised on the brane. To address this 
issue, one needs to assume fermions as propagating in the bulk (i.e. spinor fields functionally 
dependent on a) in the five dimensional warped background geometry. It has been shown 
p. Il7j| that massless fermions in a warped geometry with a RS warp factor (proportional to 
a) cannot be localised. However with additional couplings with, say, a scalar (e.g. Yukawa 
couplings) it is possible to obtain localised massless chiral fermions on the brane following 
the method developed in (in 4D flat spacetime). The question of localisation of massive 
fermions is a problem which we shall address here for the toy model with a sine-Gordon 
potential and a Yukawa coupling between the fermion and the scalar fields. 

The simplest set up for localisation is as follows. As mentioned earlier, one can achieve 
localised fermions on the brane by introducing the Yukawa coupling between fermion field 
and the scalar field, ?7i?^'F($)\l/ where F($) is a function of the rescaled scalar field (i.e. 
$ = ^). Let us first obtain the equation that governs the bulk motion of the fermions. To 
start with, we find that the Lagrangian for a Dirac fermion propagating in a five dimensional 
space with the metric (2.8), is given by 



where g = det{gab), is the determinant of full five dimensional metric. The matrices 
pa _ ^g/{cr)^/i^ — «7^) provide a four dimensional representation of the Dirac matrices in five 
dimensional curved space. Where and 7^ are the usual four dimensional Dirac matrices 
in chiral representation. The Clifford algebra in curved space, {r*^, F^} = 2g"'^ is obeyed by 
the five dimensional gamma matrices as their counterparts follow in four dimensions. These 
representations have the useful property that they can induce a chiral particle theory on the 




(3.1) 
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brane jisl . 

The covariant derivative in 5D curved space can also be calculated for the metric given 
in Eqn. (2.8) 

-D, = {d, - i/'(a)e-^('^)r^r^); V, = (3.2) 

Applying the above configuration to the Eqn. (3.1) one obtains the Dirac Lagrangian in 
5D curved spacetime in the following form 

V^Coirac = e-^^('^)^ [lef^'^h^d^ + 7^(54 - 2/'(a)) - r/^F($)] ^ (3.3) 

The dimensional reduction from 5D to 4D is performed in such a way that the standard 
four dimensional chiral particle theory is recovered. The five dimensional spinor can be 
decomposed into four dimensional and fifth dimensional parts: '^{x^, a) = '^{x'^)^{o'). Since 
the four dimensional massive fermions require both the left and right chiralities it is conve- 
nient to organise the spinors with respect to and \E'/j which represent four component 
spinors living in five dimensions given by '^l,r = |(1 -F Ts)^- Hence the full 5D spinor can 
be split in the following way 

v^(x^ a) = i^ilLixnU^) + ^nixnU^)) (3.4) 
where we have to determine ^(cr) which is an eigenfunction dependent only on a. Let us 
assume that the functions C,l.r{<^) satisfy the following eigenvalue equations 



e-^('^) [d, - 2f{a) - r]FFm U^) = -m^cr) (3.5) 
e"^^'^) [d, - 2f\a) + VFnn U^) = Ki?(^^) (3.6) 

The full 5D action then reduces to the standard four dimensional action for the massive 
chiral fermions (with mass m) , when integrated over the extra dimension Q], provided (a) 
the above equations are satisfied by the bulk fermions and (b) the following orthonormality 
conditions are obeyed. 

/oo rdo 
e-'^^^^^LMJ^ = / e-'^'^-^U^UJcr = 6mn (3.7) 
-oo J —oo 

/oo 
e-'^^'-^^L^UJcr = (3.8) 
-oo 
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The Yukawa coupling between the scalar and the fermion, with the kink solution for the 
scalar, is necessarily like an effective, variable, 5D mass for the fermions ^. This is largely 
responsible for generating the massive fermion modes in four dimensions. The dynamical 
features of the model can thus be obtained from the solutions of the eigenvalue equations 
(3.5) and (3.6). 

Let us first focus on massless (i.e. m = 0) fermions for a Yukawa coupling of the form 
?7i?^f$\l/ (i.e. F($) = $). In this case Eqn. (3.5) and (3.6) reduce into two decoupled 
equations. The solutions of which are asymptotically in the following form 



a(a) 



(3.9) 
(3.10) 



Eqn. (3.9) yields the localisation of left chiral fermions on the brane so long as, rip > 
We can as well achieve the bounded right chiral states by considering the interaction 
with an anti-kink sine-Gordon profile. 

We now turn our attention towards the massive fermions i.e. m ^ 0. Defining ^l.r{o') = 



e 2 



^l.r{<^) and making use of da^^R^a), C,r{<j) from Eqn. (3.5) and (3.6) we obtian the 



following second order, decoupled equation for the left chiral fermions whereas the right 
chiral states can be completely defined from the later equation with a prior knowledge to 
a(a) : 



dF {"(rr) 
r/^/'(cr)F($)+r/p— $'(a)-^ 



m 



Idcy) = 
a (a) 



(3.11) 
(3.12) 



We note that the equation satisfied by ^l(o") (Eqn. 3.11) is a second order differential 
equation that looks like a time-independent Schrodinger equation (more precisely like a 
'zero energy' Schrodinger equation), when recast in the following form 



m^e^/('^) + ^ + VF%<^'{a) - (,,F($) + 



Ucr) = (3.13) 
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where, |l((7) = e-f^^^iLia). 

We are interested in localisation of massive fermions around the brane. The effective 
potential for the left chiral fermions obtained from the Eqn. (3.13) should have a minimum 
at the brane in order to ensure the confinement of fermions. The Z2 symmetry of the bulk 
spacetime with respect to the position of the brane (at a = 0) puts a restriction on the 
choice of the form of the Yukawa coupling. If we demand such a symmetry for the effective 
potential -F($) should necessarily be an odd function of For example, if we consider the 
couphng as rjp^^^ (i.e. F($) = $) one obtains a symmetric effective potential given by 



Veffiy) 



27irtan"^(e'0 - 



^ + ijitanh(y) 



a 



sech(y) — //^ cosh 



(3.14) 



The above relation is written using a dimensionless variable and dimensionless parameters 



y = <^y\-^ 



m 



and for simplicity — ^ — y has been renamed as y. 

a^J\A\ 



If 



Vf 



(3.15) 




FIG. 2: For F($) = the effective potential V is plotted as a function of the rescaled extra 
dimension, y for different values of jf- The black curve corresponds to 7^ = 1, whereas the dashed 
curve represents V for 7_f = |. The parameters are chosen as = 1 and A = —6. 

Similarly, for F($) — sin$ the Z2 symmetric effective potential turns out to be : 
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FIG. 3: The effective potential V is plotted as a function of the rescaled extra dimension, y with 
different values of 7ir for F($) = sin The black curve corresponds to '^p = 1, whereas the dashed 
curve represents V for 7^ = 1/2. The parameters are chosen as = 1 and A = —6. 

As mentioned before, the crucial point here is the antisymmetric nature of the function 
F($) under the transformation $ — $. It can be be easily shown that a symmetric F($), 
such as cos$, $^ etc. yields an asymmetric effective potential. For the two different types 
of Yukawa interactions the effective potential acting on the left chiral fermions is shown 
in Fig. (3) and (4) for varying values of ^p- We have chosen a value of which may not 
correspond to the bound state. All the plots show a minimum at y = 0, which indicates 
that the massive fermions can be trapped on the brane. We notice, in addition, that the 
potential is, in general symmetric about the point = 0. Furthermore, it is worth noting 
that the value of the potential at y = in both the cases is given by : 

However, in either case, with increasing 7^, the depth of the well at the location of 
the brane becomes larger, demonstrating that more trapped states are possible for stronger 
Yukawa coupling. Additionally, the depth also increases for small v (or large p). On the 
other hand, for small 7^ the depth is less, and we expect very few trapped states (maybe 
none). The nature of the potential has prompted researchers to name it as a 'volcano' 
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potential. A fair amount of work related to bound states and resonances in such volcano 
potentials has appeared in the literature in the recent past [ist - 

We can obtain an exact solution of the equation (3.13) for 7^? = (^^ ~ exp (^^/A2r) ^j^gj^g ^ 
and A2 are defined later). However, this turns out to be an asymptotically growing function 
(recall that C, = e^^'^^'i, and we need to consider f or ^ and not ^). The function is also not 
normalisable according to the orthonormality conditions given in Eqn. (3.7). Based on all 
the above arguments we may therefore conclude that the confinement of massive fermions on 
the brane cannot be achieved without a strong Yukawa interaction. Moreover the bound on 



the coupling constant, 7]f > ^y|A|/6, required for the confinement of massless fermions also 
suggests that a large 'jp is crucial for the existence of massive bound states on the brane. 

We may also view the appearance and total number of bound states as a result of a 
competition between gravity and the Yukawa interaction. Without gravity, in a flat five 
dimensional spacetime the effective potential reduces to the Posch- Teller potential, for which 
exact bound state solutions 3] exist. But without the Yukawa interaction, we tend to get 
very few bound states and mostly quasilocalised ones which may not be normalisable. Thus, 
the presence of gravity, in some sense, destroys the bound states which were there and 
therefore, a strong Yukawa interaction is necessary in order to get them back. We shall now 
discuss these issues in some more detail. 



A. The mass spectrum of the low-lying states for F(<^) = $ and sin<^ 

We expand the effective potential about y = and retain terms upto order y^. This is the 
harmonic oscillator approximation in the neighborhood of the 3-brane. The second order 
differential equation for left chiral fermions in both the cases turns out to be : 



where we have for F($) = $: 



Uy) = 



(3.18) 



If + 




(3.19) 
(3.20) 
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^Ki 1 



k = 2\ C2- ^i—^=\ (3.21) 



and for F($) = sin $, Ci is the same as before and C2 is given by : 



From the harmonic oscillator approximation we can estimate the allowed energy levels 
and thereby find out the allowed fermion mass spectrum. The full five dimensional function 
\&(x^,cr) will have a definite parity because of the overall symmetry of the problem under 
a — > —a. Thus the right and left chiral wavefunctions will be either odd or even. We note 
here that only the even states will survive in the spectrum because of the required matching 
of the wave-function and its derivatives at the brane location (see a discussion on this in 
jl^). The possible values of /x^ are hidden in the expression : 



^^l = C, + (2n + 1) - 1^1-^=^ (3.23) 
\ a ^6|A| 

where n = 0,2,4 

The discrete mass spectrum for the trapped fermions can be obtained clearly from the 
following relations 

For, F{^) = $ : 



-7f- 



1 



(2n 



uV6 I2u 
For, F($) = sin$ : 



+ (n + 



f2n + l 



+ 



27, 



3V6Z/2 



+ 



+ 



If 



27i 

3z/2 ' 3v^, 



+ 



36z/3 



(3.24) 



/^n = -7f- 



'V6 



(2n+ 1)^ 



12z/ 



+ [n + 



1 ( 2n + l 
2^ l~6^ 



+ 



27, 



3V6Z/2 



+ 



27i 
3z/2 



+ 



27i 



3V6Z/3 



(3.25) 



where 7f = 7f + 

Notice the presence of the two parameters '-jp and v. For 7^? = we find that the higher 
states quickly approach the continuum. This is evident from the plot of /i^ versus n in Fig. 
(5) which shows that the difference between /x^ and nl^j^i tends to saturate after the first 
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0.25 

0.2- 
0.15- 

C 

0.1- 

0.05- 

- 

2 4 6 8 10 

n 

FIG. 4: The plot shows the variation of /x^ vs. n as given in Eqn. (5.34). Mass spectrum saturates 
very quickly for 7^ = 0, which has been shown for u = 1/2 (dashed curve) and 1 (continuous 
curve). The discrete levels reach to the continuum after a very few modes. 

few states. This behaviour persists for nonzero but small for which //^ goes as ^/Tf. 
However, for large 7^? the /x^ goes as '-yp and more bound states become available within the 
well. 




1. Analysis for 7^ >> 1 and F($) = $ 

We now look for bounds on the trapped masses and the number of localised states in the 
hmit 7f >> 1. The mass spectrum for this case takes the form : 

= -^J^ (3-26) 

There is an upper limit on the mass of the fermions which are trapped by the effective 
potential (3.14) because it is necessary to have a potential barrier in order to have bound 
states. Let us check the asymptotic behaviour of K//. For large values of the extra dimension 
the effective potential reduces to : 

y.„^_J^=e*_(2|^ + _i_)=) (3.27) 

To have a potential which is a continuous function of y one can match the approximate 
effective potentials given in eqn. (3.28) and (3.37) and obtain the point of continuity to be 
at y = Um ^ l^vl^ Qj^Q finds the maximum value of trapped mass from the expression 
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(3.37) using the following criterion for the vanishing of the barrier 



VeffilJ'max,ym) ~ 

Thus for >> 1, the maximum possible mass of the bound state is : 



(3.28) 



Exp 



na 



2ki 




^exp-(-^) 



(3.29) 



For example, choosing the parameters as given in Fig. (3) the value of ^JLmax ~ 0.337^. 
The maximum number of discrete massive states localised around the brane is obtained 
easily from eqn. (3.36) by introducing an extra | factor to account for the even states only. 
For a particular choice of the coupling constant and suitable range of parameters the value 
of Timax can be calculated from the integer obtained by the following relation : 



^max — lilt 



7ir7r2 aJ|A| 



I 


A 










= Int 


i 6 ' 





\/67F7r^ 



i/exp — (ttz/) 



(3.30) 



2. Analysis for » 1 and F(#) = sin$ 

Let us now look at the case with F($) = sin<l>. The mass spectrum for very strong 
Yukawa coupling reduces to the form given in Eqn. (3.36). Far away from the brane the 
Effective potential turns out to be 



2y 



eff 



- /i ev^ - {'yp + 



1 



2^6 



(3.31) 



Following similar steps, one can easily obtain the potential barrier just disappears at the 



location, y 



WW 



along the extra dimension. Applying the matching condition given in 



Eqn. (3.38) one can readily obtain /imax and nmax, which turn out to be : 



fJ'r. 



— 7i^Exp 



n^ 



Int 



' a^\M^ . 2a /|AI 
7f— j Exp W — 



a 
= Int 



6 



\/67f 



7Fexp (-i/) 



uexp (— 2i/) 



(3.32) 
(3.33) 



In Table (II) we show some values of /imax and Umax with increasing strength of Yukawa 
interaction {'jp) and varying warp-factor parameter u. 
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Yukawa 


Warping 


7f = 


= 50 


If = 


: 100 


7f = 


150 


Coupling 


Parameter, u 
















0.5 


35.81 


8 


71.62 


16 


107.43 


24 


F(^) = $ 


1 


16.33 


3 


32.65 


7 


48.98 


10 




1.25 


11.03 


2 


22.06 


4 


33.09 


6 




2 


3.39 





6.78 





10.18 


1 




0.5 


30.32 


6 


60.65 


11 


90.98 


17 




1 


18.39 


4 


36.79 


8 


55.18 


12 




1.25 


14.32 


3 


28.65 


6 


42.98 


9 




2 


6.76 


1 


13.53 


2 


20.30 


3 



TABLE I: The maximum bound state mass, fimax and number of states, rimax are tabulated using 
v = 0.5, 1 and 1.25 for increasing values of 7^ for the two different types of Yukawa coupling. 
As the coupling becomes stronger the value of iJ,max and Umax increases, which imply that the 
possibility of localisation is more pronounced for larger jp. 



Prom the above numerical values for /jimax and Umax we can make some comparative 
statements on the role of the two different Yukawa couplings. 

(i) Equating the fimax or rimax for the two Yukawa couplings we obtain v = .8. Thus for 
u < .8 the $ coupling has a chance of having more bound states compared to the sin $ one. 
The converse is true for p > .8. Note that for u = 1 (which is close to 0.8) the difference 
between the values of /imax and Umax for the $ and sin $ couplings is small. 

(ii) In the = $ case, the factor of | in the expressions for /imax and rimax and 
a corresponding factor of one for the F[$] = sin$ case seems to suggest the following 
generalisation for arbitrary : 



n^ 



Int 



= MFMoo)]} exp (-{F[$(oo)]}z/) 
({F[$(oo)]}exp(-{F[$(oo)]}i/))^ 



(3.34) 
(3.35) 



where $(00) is the asymptotic value of the kink profile, which, in our case here is 7r/2. 
Thus, in addition to the parameters 7^ and 1/, one can tune the number of bound states 



16 



using different choices for F[$] too. 

(iii) Also note that ior 1/ — 2 or larger the rimax and /imax both become relatively smaller. 
However, as we shall see later, the hfetime of these states (with large 7^ and nu > 2) is 
exponentially large. Therefore, they seem to be better candidates for describing localised 
massive fermions. 



IV. BEYOND THE H.O. APPROXIMATION, THE LIFETIME OF QUASI- 
BOUND STATES 

In the above section we dealt with the harmonic oscillator approximation to get the 
spectrum of the lowest states. The spectrum however will be corrected by the presence of 
higher order terms. In particular, as we shall sec below, the energy eigenvalues are in general 
complex with the imaginary part being related to the lifetime of the quasibound state. 

To analyse this in some generahty let us go back to the original Schrodinger like equation 
written using the variable a. One may easily recast things using y or y. Expanding all terms 
in Taylor series upto order cr^ we obtain 

dULia) + [m^ + g'\,=o - {aia' - a2a^}] ^(a) = (4.1) 

where. 



02 



2 HI 



m 



9° + ^ 



J (7=0 



' j-//// y//2 ^ 



+ 



9" 



g'g'" 



4! 



(4.2) 
(4.3) 



a=0 



In the above expressions the prime denotes derivative with respect to a and g{a) — 
77fF($) +f'{a) /2. For the constant ai and 02 positive we can treat the potential as harmonic 
oscillator potential plus a perturbation term proportional to if the parameter 02 can be 
assumed as very small compared to ai. Given the symmetry of the system under the 
transformation a —a wc put a restriction on the Yukawa coupling term F{^) such that 
the even order derivatives of g{a) vanish at the location of the brane. In addition, keeping 
in mind the problem in hand, the odd order derivatives of f((7) is also assumed to be zero 
at the brane location. It is worth noting that under these assumptions and within the 
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approximation the potential has a volcano profile. Keeping further terms will not alter the 
nature of the potential to a very great extent. Let us now find out the lifetime of the quasi- 
bound state for the particular choice of coupling = sin$. All the above mentioned 
properties are satisfied by the specific f(cr) and g{a) discussed earlier in this article. We 
assume |A| = 6 for convenience. In particular, we have, 

g{a) = (^rjp + ^) tanh (^^) (4.4) 
Using the general formulae for ai and 02 we obtain: 



ai 



rri^ ^ rfp ^ r]F 
1/ 1/2 z/^ 



a2 



m 



2i/2 



V 



(4.5) 
(4.6) 



"iv) ' 3 z/4 

where (vf + is replaced hy fjF- It is useful to mention here that factors of sqrt\A\/6 
which have been taken to be one with the choice |A| = 6 provide the correct dimensions for 
each term in the above and in the subsequent discussions. We note that in the range of the 
warping parameter u > 2 the perturbation constant 02 is reasonably small compared to a. 
One can see this by further approximating the potential for moderately large u and large 
f]F- Ignoring terms proportional to and f]F we find that the effective potential can be 
written as : 



V{a) 



Using the above approximate potential we can use the method developed in [IJj for the 



(4.7) 



calculation of the lifetime, r 



1 



2ImE 



(where ImE is a small negative quantity). We briefly 



summarize their method below. 

The imaginary part ImE is given by the formula : 



ImE 



d(T2 



(4.8) 



The denominator is calculated using the zeroth order wavefunction (harmonic oscillator) 
in Rayleigh Schrodinger perturbation theory. In such an approximation we assume that the 
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dominant contribution to the norm comes from the small values of a. For the numerator, we 
use the WKB method with the assumption that the probability current is dominated by the 
contribution of the wave function in the classically forbidden region. Same normalisation 
for the small and large a solution has been guaranteed by the asymptotic matching of the 
respective solutions in these regions. 



It is easy to show (following 



16|) that if riF is large and u >2 the lifetime of the metastable 



lowest state of the harmonic oscillator goes as 

1 

CT — exp(r/Fi/) (4.9) 
Vf 

From the above formula it is clear that the lifetime is exponentially enhanced for large 
fipu. Note that i/ > 2 is required in order to treat the term as a perturbation. If f]pv>bQ 
then the lifetime will be close to the age of the universe (lO^Oyears). Therefore large f]F 
would mean values greater than 25 (for v = 2). For such values of the parameters the 
massive mode can be, more or less assumed to be stable. 

One can also evaluate the lifetime for the $ Yukawa coupling though qualitatively the 
results will not differ except for factors here and there. 



V. CONCLUSIONS 



Let us now summarise pointwise the results obtained in this article and briefly discuss 
some possible avenues of future research. 

(i) In the presence of a sine-Gordon potential we have first discussed an exact solution 
for the warp factor and the scalar kink by solving the full set of Einstein-scalar equations 
in the presence of a negative cosmological constant. The background spacetime however, 
in this case, has a non-constant Ricci scalar which is asymptotically negative. The matter 
stress energy which generates this warped line element satisfies the null energy condition 
with isotropic negative pressure in the vicinity of the brane. The pressure along the extra 
dimension is, however, zero at the brane location and positive away from the brane. 

(ii) Fermion fields have been studied in the background geometry of the exact solution 
mentioned in (i) above. We couple the scalar field to the spinor field through a Yukawa 
coupling and write down, following standard methods, the effective Schrodinger like equation 
for massive fermions. Two different types of Yukawa interactions are considered here, for 
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each of which we analyse the issue of locahsation of massive fermions. The difference between 
the effective potentials of the two Yukawa interactions lie in the shape and depth of the well 
at the brane location — for the same value of the Yukawa coupling parameter the one 
which is more steep is prone to have more bound states. Using approximate methods we 
are able to obtain the maximum mass, the low-lying mass spectrum and the maximum 
number of trapped states for each case, for large values of the Yukawa coupling parameter. 
We also calculate the lifetime of these quasilocalised states and show that it can indeed be 
exponentially large for large 7^ and moderate values of v. In addition, we also show that 
for zero or small values of 7^ (fixed p) very few bound states may exist and the spectrum 
becomes almost continuous beyond these discrete states — a fact demonstrated through the 
vanishing of the level spacing for large values of the quantum number n. 

For massless fermions, we verify previously known results. Only the left or the right chiral 
fermions can be localised under special restrictions on the value of the Yukawa coupling 
constant and other parameters that appear in the warp factor/kink solution. 

(iii) An obvious question related to the background geometry arrived at in this paper is : 
is gravity localised on the brane in this model? The answer to this question, as far as we have 
been able to analyse, is in the affirmative. The potential felt by the graviton fluctuations is 
similar to the one in the RSII model with a well in the neighborhood of the brane followed 
by a barrier and a decay to zero as we go deep into the bulk. We shall elaborate on the 
localisation of gravity in our model in a future publication. 

Acknowledgment s 

SK thanks S. Randjbar Daemi for discussions and the ASICTP, Trieste, Italy for support 
during part of the period when this work was done. RK thanks IIT Kharagpur for support. 



[1] T. Kaluza, Sitzungsber. Preuss. Akad. Wiss., Phys-Math.Kl., Berlin Math. Phys., Bd. Kl 966 

(1921) ; O. Klein, Z. Phys. 37, 895 (1926) 
[2] A. Salam and J. Strathdee, Ann. Phys. 141, 316 (1982); M. J. Duff, B.E.W. Nilsson and C. N. 

Pope, Phys. Rep. 130 C, 1 (1986); T. Appelquist, A. Chodos and P. G. O. Preund, Modern 



20 



Kaluza-Klein Theories, Reading, MA, Addison- Wesley (1987) 
[3] M. S. Green, J. H. Schwarz and E. Witten, Superstring theory, Cambridge University Press 

(1987); J. Polchinsky, String theory, Cambridge University Press (1997); 
[4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999); ibidPhys. Rev. Lett. 83, 4690 

(1999) 

[5] Y. Kim, C. O. Lee, I. Lee, J. J. Lee, |hep-th/0307"023l A. Padilla, |hep-th/0210'217l R. Dick, 
Class. Quant. Grav. 18, Rl (2001); C. Csaki, |hep^/0404096l 

[6] B. Bajc and G. Gabadadze, Phys. Lett. B 474, 282 (2000); S. L. Dubovsky, V. A. Rubakov 
and P. G. Tinyakov, Phys. Rev. D 62, 105011 (2000) 

[7] M. Besancon, Experimental introduction to extra dimensions, hep-th/Q106165| ; Y. A. 
Kubyshin, Models with extra dimensions and their phenomenology hep-ph/0111027| and refer- 
ences therein. 

[8] D. Langlois, ^gr-qc/0102007 (Proceedings of the 9th Marcel Grossmann meeting, July 2000 
(Rome); gr-qc/0205004| (Proceedings of Journees Relativistes, Dublin 2001 and references 
therein ; E. Flannagan, S.-H. Henry Tye and I. Wasserman, Phys. Rev. D 62, 024011 (2000); 
H. Stoica, S.-H. Henry Tye and I. Wasserman, Phys. Letts. B 482, 205 (2000); J. S. Alcaniz 
, astro- ph/0202492l 

[9] W. D. Goldberger and M. B. Wise, Phys. Rev. D 60, 107505 (1999); H. Davoudiasl, J.L. 
Hewett and T.G. Rizzo, Phys.Lett. B473 43 (2000); S. C. Davis, JHEP 0203, 054 (2002); E. 
E. Flanagan, S.-H. Henry Tye and I. Wasserman, Phys.Lett. B522 155 (2001); B. Mukhopad- 
hyaya, S. Sen, S. SenGupta, Phys.Rev. D65 124021 (2002); M. Gogberashvih, Phys.Lett. 
B553, 284 (2003); Rui Neves, Cenalo Vaz, hep-th/0302030 
[10] C. Ringeval, P. Peter, J. P. Uzan, Phys. Rev. D 65, 044416 (2002); S. Ichinose, Phys.Rev. 
D66, 104015 (2002) 

[11] For a useful summary on the details of the energy conditions see Lorentzian Wormholes : from 
Einstein to Hawking by Matt Visser (AIP, 1995) 

[12] M. Gremm, Phys. Lett. B478, 434 (2000); Phys. Rev. D62,044017 (2000); K. Ghoroku and 
M. Yahiro, |hep-th/0305150t A. Kehagias and K. Tamvakis, Mod. Phys. Lett. A 17, 1737 
{2002),ibid. Phys, Letts. B504, 38(2001); M. Giovannini, Phys. Rev. D 64, 0106041 (2001), 
Phys. Rev. D 65, 064008 (2002); S. Kobayashi, K. Koyama and J. Soda, Phys. Rev. D65, 
064014 (2002) 



21 



[13] O. DeWolfe, D. Z. Friedman. S. S. Gubser and A. Karch, Phys. Rev. D62 046008 (2000); 

C. Csaki, J. Ehrlich, T.J. Hollowood and Y. Shirman, Nucl. Phys. B581, 309 (2000) ; S. 

Ichinose, Phys. Rev. D 65,084038 (2002) ibid. Class.Quant.Grav. 18, 5239 (2001); A. Campos, 

Phys.Rev.Lett. 88 141 (2002) 
[14] M. Cvetic and H. Soleng, Phys. Rept. 282, 159 (1997); P.M. Cowdah, H. Lu, C.N. Pope, K.S. 

Stehe, P.K. Townsend, Nucl.Phys. B486 49 (1997); K. Skenderis and P. Townsend, Phys. 

Letts. B468, 46 (1999); K. Behrndt, Nucl.Phys. B573 127 (2000); i K. Behrndt, M. Cvetic, 

Phys.Lett. B475 253 (2000); M. Cvetic, H. Lu, C.N. Pope, Class.Quant.Grav. 17 4867 (2000); 

C. Csaki, J. Erlich, C. Grojcan and T. Hollowood, Nucl. Phys. B584, 359 (2000); K. Behrndt, 

C. Herrmann, J. Louis, S. Thomas, JHEP 0101 Oil (2001); G. L. Cardoso, G. Dah'Agata, D. 

Lust, JHEP 0107 026 (2001); A. H. Chamseddine, W. A. Sabra Nucl.Phys. 630 326 (2002); 

CM. Hull, JHEP 0111 061 (2001); M. Bianchi, D. Z. Preedman, K. Skenderis , Nucl.Phys. 

B631, 159 (2002); K. Behrndt, M. Cvetic, Phys.Rev. D65 126007 (2002); L. Anguelova, C.L 

Lazaroiu, JHEP 0209 053 (2002); S. L. Cacciatori, D. Klemm, W. A. Sabra, JHEP 0303 023 

(2003) 

[15] E Caliceti, V. Grecchi and M. Maioh, Commun. Math. Phys. 176,1 (1996), J. Zamastil, J. 

Cizek and L. Scala, Phys. Rev. Letts. 84, 5683 (2000); F. J. Gomez and L. Sesma, Phys. Letts. 

A 301, 184 (2002);; M. M. Nicto, Phys. Letts. B486,414(2000) 
[16] C. M. Bender and T. S. Wu, Phys. Rev. 184, 5 (1969); C. M. Bender and T. S. Wu, Phys. 

Rev. D 7, 6 (1973); J. Zamastil, V. Spirko, J. Cizek, L. Skala and O. Bludsky, Phys. Rev. A 

64, 042101 (2001) 

[17] S. Randjbar-Daemi and M. Shaposhnikov, Phys. Letts. B492, 361-364 (2000) 
[18] R. Jackiw and C. Rebbi, Phys. Rev. D13,3398 (1976); Y. Grossman, N. Neubert, Phys. Lett. 
B474 361 (2000). 

[19] L.D. Landau and E.M. Lifshitz, Quantum Mechanics, Course of Theoretical Physics, Vol. 3, 

Third Edition (Butterworth bitemwb 
[20] S. Weinberg, Gravitation and Cosmology (John Wiley and Sons, 1971) 



22 



